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Abstract—Parametrical study for performances comparison between classic orifice pulse tube and hybrid Stirling pulse 
tube refrigerators. An internal dynamic comparative study of orifice pulse tube (OPT) and Stirling hybrid pulse tube (HPT) is 
made, by thermodynamic modelling of the two systems. Basical structure of the model is established based on mass and energy 
conservation balances, written for separate control volumes of complete refrigerator. This method makes clear typical common 
parameters of the two mechanisms and specific importance of pistons phase shift and swept-volumes ratio for HPT. Starting from 
identical numeric characteristics for the two refrigerators, the study considers principally phases' evolution for pressure and mass 
flow rate at tube extremities. Greatest setting easiness, lower dimension for equivalent performances, give to the Stirling hybrid 
pulse tube good chances to be developed in miniature cryocooler area. © Elsevier, Paris. 

Stirling cryocooler / pulse tube / comparison / phase shift / parametrical study 

Resume — Une etude comparative de la dynamique interne des refroidisseurs de type pulse tube a reservoir et orifice (OPT) et 
Stirling-pulse tube (hybride : HPT) est menee de maniere theorique grace a la modelisation thermodynamique de chacun des deux 
systemes. Le schema de base du modele est etabli a partir des bilans de masse et d'energie ecrits pour les differents volumes de 
controle du refroidisseur complet. Cette faijon de proceder fait ressortir les parametres communs au fonctionnement des deux 
appareils ; le role du dephasage des pistons et du ratio des cylindrees est mis en evidence pour le HPT. A partir d’un cas de figure 
numerique avec des donnees communes, I'etude est surtout axee sur les phases entre la pression interne et les debits massiques 
aux extremites du tube. Une plus grande souplesse de reglage et une taille moindre pour des performances equivalentes font du 
HPT un systeme apte a se developper dans le domaine des microrefroidisseurs. © Elsevier. Paris. 
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Nomenclature 

c v ,c p capacite calorifique isochore et 


• r x * _ 1_1 

isobare. J-kg 1 K -1 

d diametre capillaire. m 

E efficacite du regenerateur 

Eq parametre modele thermique 

ref. [15] 

h enthalpie massique. J-kg -1 

H flux d 'enthalpie. W 

Ko= facteur de perte de charge. kg-s -1 'Pa -1 

La longueur du regenerateur. m 

l longueur capillaire. m 

m masse de gaz. kg 


rri 

debit massique de gaz. 

kg-s- 1 

p 

pression inoyenne du gaz. 

uj d 2 

Pa 

Pe = 

-nombre de Peelet 

4 a 


Q 

quantite de chaleur. 

J 

r 

constante des gaz parfaits. 

J-kg-'-K -1 

s 

section. 

m 2 

t 

temps. 

s 

T 

temperature. 

K 

u 

vitessc . 

m-s -1 

u 

energie interne. 

J 

Vo 

volume mort du cylindre au PMH 

m 3 

w 

travail mecanique des forces de 



pression.. 

J 

X 

variable d'espace. 

m 
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Indices 


a 

ambiante 

c 

cote chaud du tube 

cp 

compresseur 

d 

regenerateur 

di 

en entree 

d 2 

cn sortie 

e 

en entree ou cote froid 

ec 

echangeur de refroidissement 

ge 

gaz entree du regenerateur 

gs 

gaz sortie du regenerateur 

rep 

dephaseur -det endear 

s 

en sortie 

se 

zone froide du tube 

sc 

cylindre de compresseur 

t 

Expos ants 

tube 


* irulique la dependance avec le 

signe du debit gazeux 

Lettres grecques 

3 phase models themiique ref. [15] 

c p 

7 = — coefficient isentropique du gaz 
Cy 

ft viscosite dynatnique. 

p masse volumique. 

ui pulsation. 


kg-m 1 -s 1 
kg-m'"* 


1. PRINCIPE DU REFRIGERATEUR 
HYBRIDE STIRLING PULSE TUBE 


Un refrigerateur a tube a gaz pulse (figure 1 ) 
coinporte generalement un volume de compression defini 
par un piston mobile dans son cylindre. un echangeur 
de refroidissement, un regenerateur separe du tube a 
gaz par la zone froide, 1 c tube lui-mcme, un echangeur 
chaud et un orifice calibre mettant en communication 
Lextremite chaude refroidie du tube avec un reservoir 
suffisainment grand. Le systeme dit «a deux orifices». 
coinporte line derivation partielle du gaz entre l’entree 
du regenerateur et rextremite chaude du tube et 
se montre plus efficace au niveau de la puissance 
frigorifique et des temperatures minimales obtenues, car 
la reduction du debit le traversant ameliorc Lefficacite 
du regenerateur. 

Pour comprendre la dynamique interne du systeme 
du refrigerateur a tube a gaz pulse, on peut en resuiner 
le fonctionnement de la maniere suivante. 

Le gaz. comprime lorsquc le piston se deplace pour 
reduire le volume de la chambre, subit unc compression 
adiabatique : il est ramene a. la temperature ambiante 
grace a Lechangeur refroidisseur dispose en sortie. Une 
partie du fluide traverse le regenerateur pour entrer 
dans le tube, une autre partie sort du tube pour 



Figure 1. Diagramme schematique devolution du flux d’en- 
thalpie dans les refroidisseurs Stirling, pulse tube a orifice et 
Stirling pulse tube hybride. 

Figure 1. Simplified diagram of enthalpic flow into Stirling 
orifice pulse tube and hybrid pulse tube refrigerators. 

penctrer dans le reservoir. Un echangeur, situe entre 
rextremite du tube et le reservoir, permet de ramener 
ce gaz a la temperature ambiante et evacue une energie 
correspondant pratiquement a la chaleur soutiree par 
Lechangeur froid. Lorsque le piston repart en sens 
inverse, la detente refroidit le gaz dans le tube a une 
temperature inferieure a celle qu'il avait en sortant du 
regenerateur. puis la pression diminuant encore, elle 
devient inferieure a celle du reservoir : le gaz qui on 
ressort repousse alors le gaz froid, qui retraverse la zone 
a refroidir et le regenerateur en sens inverse. 

On explique Lexistence d'une production frigorifique 
a Lextremite « froide» du tube par la creation d’un flux 
(Lenthalpie axial qui traverse celui-ci [1-3]. Ce flux, dont 
la moyenne temporelle s’exprime par : 

, 

(H) = ^ J J pu(t)hdS 

r 2 ? r 

-sfhii //(')“««' 15 <‘» 

est. maximal dans le cas de variations periodiques de 
la pression et de la vitesse du gaz. lorscjue ces deux 
grandeurs sont en phase. On doit done realiser un 
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dephasage adequat de (Jt/2) entre la pression du gaz et 
son deplacement dans le tube (integration de u(t)), pour 
obtenir une production frigorifique inaxirnale. Dans une 
machine Stirling, ce dephasage est, regie par le piston 
deplaceur : dans le tube a orifice, il est regie par 1 a. 
difference de pression avec le reservoir, et dans un tube 
ferine ou un tube tliermoacoustique par le transfert de 
chaleur irreversible avec les parois. 

La figure 1 donnc une idee de la repartition du flux 
d'enthalpie le long de ces divers refroidisseurs. On notera 
que, pour le regime periodique etabli, le gaz avant une 
vitesse rnoyenne nulle dans le regenerateur. il v subit 
une transformation quasi isotherme, car sa capacite 
thermique est negligeable vis-a-vis de celle de la matricc 
solide ; le flux d'enthalpie est alors mil. 


1 


pu(t) di = 0 


o 




cp 

tS 




pu(t)dt dS 0 (2) 


On en deduit, pour un tube a gaz pulse, les egalites : 
\W\~m et \Q' t + Qi\ = \Qc + W'\ (3) 


Les tubes a gaz pulse classiques avec orifice et 
reservoir (OPT) ajoutent a l’encombrement important 
de ce dernier le desavantage dun reglage de phase 
restreint et peu commode. Il est en effet impossible 
de vraiment agir efficaccment stir les phases, entre la 
pression dans le tube et les debits gazeux aux extremit.es. 
sans modifier egalement leurs amplitudes. Ce faisant. on 
risque de diminuer largement la dynamique interne de 
pression dans le tube, facteur crucial pour robtention 
d'une production frigorifique sufflsante. 

Le handicap du refroidisseur Stirling est double lui 
aussi : fonctionnement du piston detendeur en zone 
froide. avec tous les problemes de frottements et de 
jeux que cela impose, et niveau vibratoire important, 
occasionne par la contiguite de la zone froide avec le 
cylindre de detente. 

Le tube a gaz pulse hybride (HPT) (figure 2b) resout 
en partie ces problemes : meilleure maitrise du reglage 
de phase entre les deux pistons et report du piston 
dephaseur detendeur dans une zone a temperature 
ambiante. L’idee originale en est due a Kasuya et al. 
[4-9]. mais d'autres auteurs [10 12 ] ont depuis demontre 
experimentalement l’efficacite de ce systeme, capable 
d’atteindre actuellemeut une temperature de 23 K [13] 
avec une machine a un soul etage. 

Nous presentons une etude theorique comparative 
entre les refroidisseurs OPT et HPT. axec stir les 
differences obtenues entre les evolutions temporelles des 
pressions internes et des debits massiques aux extremities 
pour les deux appareils. Les parametres importants an 
niveau de la caracterisation des performances du HPT 
sont mis en lumiere et leur influence est discutee. 


Te me>0 


m (t) P(t) 
V(t) 


Qf 
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Ts ms<0 


echangeur 


reglage 
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Vt 
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Figure 2. a. Volume deformable utilise pour les bilans et 
notations, b. Volumes de controle definis pour le modele et 
notations utilisees. 

Figure 2. a. Deformable volume used for balances and nota¬ 
tions. b. Control volumes and notations used for modelling. 


2. SCHEMA DIRECTIF 

POUR LA MODELISATION 
PARAMETRIQUE DU HPT 


Pour la modelisation du HPT. la methode choisie 
est de caractere thermodynamique ct dccoule de celle 
que nous avons etablie precedemmcnt pour les machines 
Stirling puis pour les Pulse Tubes a deux orifices [14]. Les 
equations theoriques de base sont obtenues a partir des 
Ijilans de matiere et d'energic etablis pour un volume 
de controle ouvert de forme quelconque. La figure 2a 
illustre un volume dc base, sur la figure 2b sont designes 
les sept volumes definis separement. pour lc refroidisseur 
complct HPT : 

bilan de conservation de masse : 

m„ + r'rn = Am/At (4) 

bilan de conservation de Fenergie : 

— P d V + r T„ dm c +rT s dm s + 8Q + Cv - T e dm e 

+ cv T s dm s = c r A(tnT) 


avec la loi des gaz parfait.s : 


P V = m r T 



+ c p T e m.e + Cp T r rh s + 8Q 


d(mT) _ 1 d(PV r ) 

dt 7 — 1 df 


( 6 ) 
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Les mouvements des deux pistons sont sinusoidaux 
et presentent un dephasage reglable «phas» : 


Vcp(t) = ~ (1 - sin uit) + Vo; 

V rcp (t) = - 7 p(l - sin (cat + phas)) + Vo (7) 

v:,c et V’ rcp sont respectivement les volumes balayes des 
deux cylindres. 

Pour le regenerateur on fait l'liypothese dc variations 
lineaires de la pression et de la temperature de la matrice 
solide et du gaz a 1'interieur, on tient aussi conipte de 
son efficacite E pour le calcul de la temperature du 
fluide qui en ressort. 


P d (x,t) = Pc(t) + P c (t) 

r mat rice(xi) = ^ T * X + T & 


(8) 


soit, pour le gaz contenu dans le regenerateur : 

Td(x,t) - £r mtt trice(x,i) + Tei«rfe(l ~ E) 

= E Tc ~ T& X + E P a + Xe„t r ee(l ~ E) (9) 


La definition theorique de E impose normalement sa 
dependance avec le temps E(f) : 


(Tgs(t) 


r ge ) dt 


E(t) = 


t (Trnatrice Pge) 


Cette fonction etant pratiquement impossible a 
definir. nous sommes malheureuseinent contraints de 
faire l’hypothese restrictive d'une efficacite constante. 

Le debit massique de fluide traversant le regenerateur 
est represente par la forme suivante : 


m dl = K 0 (P c (t) - P t (t)) (10) 


la reference [15]. nous avons presente des expressions 
tie calcul ties flux thermiques periodiques sous forme 
de nombres complexes, afin de tenir compte des phases 
differentes ent.re pression. temperatures du gaz et de la 
paroi et flux parietal. La partie reelle du flux thermique 
sera ici representee par la relation : 

S t f Q, : dx = - (P r - P) cos 3 (12) 

./o vPe 7-1 

Pour la representation des debits aux jonctions on 
utilise la loi de Poiseuille pour un capillaire. donnee en 
fonction de la difference des pressions. de la longueur 
et du diametre. L'organe est caracterise par un terme 
d’ouverture d‘expression generale : 



soit m r = Kx{P rep(t) - Pt(t)) (14) 

Le tableau I resume toutes les variables du systeme 
complet d'equations obtenu. L'ensemble des bilans 
fournit un systeme tie sept equations differentielles 
principales pour les variables, 

PAt).P(t). P rrp (t), Tc„(t), T t (0,f), Tt(L.t), T lcp (t) 

A partir de ces valeurs on dispose aussi tie sept 
equations pour les calculs ties debits 

the, disc- rHcp, m t u, m r , mi. th ec 

et ties equations pour le calcul de benergie 8Qt soutiree 
a la source froitle et du travail ecliange par le gaz avec 
chaque piston 

81 L = + mcpT’ 

La resolution du systeme doit etre faite etape par 
etape en fonction des directions des differents debits 
(positifs en entree et negatifs s'ils sortent) et des 
temperatures ci-dessous. Le tableau II illustre les 8 
phases de calcul successive durant un cycle complet 
du refroidisseur. 


Pour le tube a gaz. nous supposons aussi que la 
repartition des temperatures est lineaire : 

(T* - T*) 

T t (x,t) = 1 * ■■.. . ■ & x + T; e (11) 

Dans cette relation, les valeurs des temperatures sont 
a determiner suivant les signes des debits gazeux aux 
deux extremites du tube ; on se reportera ci dessous aux 
relations (15) pour le detail. 

Le tube a gaz est le seul volume considere comme 
non isotherme ; durant les periodes de compression et de 
detente les temperatures aux extremites varient, ellcs 
aussi, de fagon periodique. Des echanges de chaleur 
alternes surviennent ent.re le gaz ct la paroi. Bien que 
localises surtout, dans la couche limite et present ant une 
moyenne temporelle nulle (stockage- -destockage). ils ne 
sont pas sans influence, notamment sur les phases. Dans 


T* — T t (L t ,t) 

si rhe > 0 

P s * e =P(o,t) 

si m so > 0 

II 

si rh, c < 0 

= Pe 

Si 777.se < 0 

Pep — Tcp(t) 

si rhep < 0 

Prep = Prep (t) 

si m r < 0 

= T a si ih c p >0 = Pa si m r > 0 

(15) 

Apres adimensionnalisation. la resolution des equa¬ 
tions est effectuee numeriquement a l’aide du logiciel 
Mathematica. 

t' = U)t 

("0 

II 

73| n "0 

P 

P' — _1 p' 

~ p r rcp 

_ Prep 

P 

, rt 1 « P 

Vcp ~ i: 

^ r t V rC p 

rcp “* Vt 

v' = — v\ = 

v ec r r v d 

Vt 

V A 
' V t 

v; 

v' - — 

r se r 7 

n 
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TABLEAU 1 / TABLE 1 

Definitions des notations des pressions, volumes et debits massiques 
pour les 7 volumes de controle utilises pour les equations de bilan. 

Definitions of notations for pressures, volumes and mass flows 
of the 7 control volumes used in balance equations. 


Compresseur 

Echangeur 

Regenerateur 

Partie froide 

Tube a gaz 

Partie chaude 

Dephaseur 

Parametre 

Vcp(f) 

V et: 

Vd 

V's, 

V, 

Vc 

Vrcp(t) 

Temperature 

Tc(t) 


Td(x,t) 

T, 

T t (x,t) 

T n 

Tcp(f) 

Masse de gaz 

m cp (f) 

Tries (t) 

■m d (t) 

m c (t) 

m t {t) 

m c (t ) 

m r (t) 

Pressions 

Pit) 

Pit) 

Pd{x,t) 

Pit) 

Pit) 

Pit) 

Prep (t) 

Debit entree 

rhcp(f) 

m cp (t) 

rhdi(t) = rh tx (t) - m\{t) 

rh d2 (t) 

m S c(f) 

ihc(t) 

m r (t) 

Debit sortie 

rh C p(f) 

m ec (f) 

rhd 2 (t) 

rh se {t) 

m c (f) 

mi{t) — m r (t) 

m r (t) 

Energie 

VVcp 

Qec. 


Q 

Q 

Qc 

tVrcp 


TABLEAU II / TABLE II 

Decoupage des phases durant un cycle complet pour I’ecriture des equations de bilan. 

Cutting of phases, during a complete cycle, used in balance equations. 

Phase 1 

Phase 2 

Phase 3 

Phase -1 

Phase 5 

Phase G 

Phase 7 

Phase 8 

Phase lb 

m cp < 0 

nicp < 0 

rricp < 0 

rhep > 0 

rricp > 0 

rricp > 0 

m cp > 0 

m cp ^ 0 

rh (: p < 0 

rile < 0 

rhe > 0 

rile > 0 

m c > 0 

rii <; > 0 

rh e < 0 

rhe < 0 

rile < 0 

rite < 0 

ritse < 0 

tYi se 0 

O 

V 

r 

ritse < 0 

ritse > 0 

ritse > 0 

riise > 0 

rhso > 0 

m s «. < 0 

rh, < 0 

m r < 0 

m r > 0 

m r > 0 

rilr > 0 

riir < 0 

tit, < 0 

m r < 0 

rri r < 0 


Cette fa^oti de proceder met en lumiere des groupc- 
ments importants dans la modelisation : il y a lieu de 
distinguer les facteurs d'ordre geometrique on thermi- 
que de ceux lies aux ecoulements au travers des orifices 
divers : 

- facteur de production frigorifique 


_ 8 Q f 

PiO V't 


(16) 


- facteurs tliermiques 
7 coefficient isentropique du gaz 

cos 3 facteur de transferts tliermiques aux parois 

JPi 

E efficacite du regenerates 

T 

~ rapport, des temperatures 

C 

( 20 ) 


- facteurs geometriques : rapports de volumes : 

Vcp Id ’ ' ' i. Vec * , 

—■- compresseur — regenerateur - echangeur 

Vt Vt Vt 

^ 7 - zone chaude rpr zone froide ^ '' r p dephaseur 

Vt 


Vt Vt 

- facteurs d'ecoulement du gaz : 

r T a K r r T a n 


AVr = 


■Vt 


ft ( Ci p \ 

pa 1 — 1 

r P V 2 / 


(17) 


(18) 


ll. Vt 8 It fj. 

facteur d orifice du reservoir 


A' k „ = 


r T a A'o r T a m d 


(19) 


; VV w V t A P 
facteur de perte de charge du regenerateur 


3. COMPARAISONS DES RESULTATS 
DE L’ETUDE PARAMETRIQUE 
DES REFRIGERATEURS A TUBE 
A GAZ CLASSIQUES ET HYBRIDES 


Les graphiques presentes dans cette etude montrent 
des evolutions temporelles parfois tres differentes d’une 
allure purement sinusoi'dale. Notons, pour la suite 
de 1 ’expose, que deux grandeurs seront par abus de 
laugage dites en phase lorsque leurs valeurs minimales 
se situeront au mane instant (pression minimale 
et temperature minimale par exemple) : la reference 
origine pour les phases est prise nulle lorsque le piston du 
compresseur est en position mediane dans sa chainbre : 


le volume de gaz est alors + Vo. 


Dans ce paragraplie. les valeurs de base choisies 
coniine exemple pour les premiers calculs concernent un 
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tube a gaz pulse a orifices. OPT, contenant. de Lhelium 
(7 = 1 , 66 ) et. sont les suivantes : 


Psc 

¥ 


= 1 , 8 ; 



V t 


14 

V, 



- pour le regenerateur : A k „ = 15 et E = 1 (valeur 
maximale) 

- temperatures des echangeurs des cotes chaud et 
froid : T c = T a = 300 K et T e = 240 K. 

- transferts theriniques alternes : 

En 

— 7 = cos 3 = 0,0141 
sfPe 

Pour comparaisons, les calculs pour le refrigerateur 
hybride HPT sont effectues avec les memes valeurs. 

Les figures 3a-3c representent les variations en 
fonction du parametre d’ouverture K kr . de la pression du 
gaz et des debits massiques d'extremite pour un OPT. 
Augmenter A kr fait diminuer la production frigorifique 
du tube, surtout parce que le rapport des pressions 

extremes max diminue. La phase de la pression diminue 

1 min 

par la meine occasion, et les deux debits s'accroissent : 



Figure 3. AT r — 0,1 ; 0,5 ; 1 ; 1,2. a. Evolution temporelle de 
la pression dans le tube. b. Evolution temporelle des debits 
massiques aux extremites du tube. c. Evolution temporelle 
du facteur de production frigorifique du tube. d. Facteur de 
production frigorifique et mecanique du tube en fonction de 

I’ouverture A" k r, cas — =0,8. 

Tc 

Figure 3. A' kr = 0.1; 0.5; 1; 1.2. a. Pressure temporal evo¬ 
lution in an orifice and reservoir pulse tube. b. Mass flow 
temporal evolution on the extremities, c. Refrigerating factor 
temporal evolution, d. Refrigerating and mechanical factors 
for the tube depending on reservoir orifice aperture AT r , 
Tf 

when — = 0.8. 

T c 


en fait, leurs phases sont assez peu affectees. On verifie 
que la configuration la plus favorable pour la production 
frigorifique est. obtenue quand la valeur de A' kr est 
petite ; pression et debits sont alors a peu pres en phase. 
Pour les valeurs elevees de A' kr . le tube peut meme 
inverser ses echanges de chaleur aux deux extremites et 
se transformer en generateur d’energie acoustique. 

Pour le HPT, deux parainetres jouent differemment 
un role actif sur les amplitudes et les phases, ce sont 

y 

le rapport des cylindrees compresseur-dephaseur —^ 

V'cp 

et. la valeur dtt dephasage phas entre pistons. Les fi¬ 
gures 4 et 5 montrent riiifluence du ratio de cylindree 
pour des valeurs respectives de phases entre pistons de 
0,125 7 t/2 et 0,5 w/2. La combinaison des deux effets est 
assez complexe. Pour les faibles dephasages, le rapport 

- p niax augmente avec le ratio ; il a tendance a decroitre 

■* min 

au-dela dame certaine valeur de phas. Lorsque Ton aug¬ 
mente le dephasage phas, on diminue aussi la phase 
entre la pression et le debit cote chaud et Ton augmente 

Lamplitude de celui-ci ; quand on augmente la 

V''cp 

phase du debit cote froid est inchangee et son amplit ude 
diminue dans le cas phas = 0.125 Jt/2. La phase diminue 
avec line amplitude rnoyenne a peu pres constante pour 
le cas phas = 0,5 7 t/ 2 . La production frigorifique sui- 
vra. un comportement semblable a celui du rapport des 




Figure 4. Ratios de cylindree 0,25; 0,5 ; 0.75; phase 
0,125 71/2. a. Evolution temporelle de la pression dans le 
tube a gaz. b. Evolution temporelle du facteur de production 
frigorifique. c. Evolution temporelle du debit massique de 
gaz a I'extremite chaude. d. Evolution temporelle du debit 
massique de gaz a I’extremite froide. 

Figure 4. Swept-volume ratios 0.25, 0.5, 0.75 ; phase shift 
0.125 71/2. a. Pressure temporal evolution inside tube of HPT 
refrigerator, c. Mass flow temporal evolution at hot extremity 
of HPT refrigerator, d. Mass flow temporal evolution at cold 
extremity of HPT refrigerator. 
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pressions : pour chaque ratio de cylindree il exist.e une 
valeur optimale a choisir pour le dephasage (figure 6). 

La figure 7 illustre 1 : influence de phas pour un ratio 
de cylindree valant 0,25 ; on notera la possibilite de 
regler le debit cote chaud sans trop influencer celui du 
cote froid. ce qui n'est pas possible avec le tube a orifice. 


4. CONCLUSIONS 

En raison des hypotheses simplificatrices adoptees, 
l'objectif de comparaison des deux types de refrigera¬ 
teurs est atteint sur la base d'une modelisation thermo- 
dynamique de type “premier ordre». 




Figure 7. Phases 0,125 7 t /2 a 0,8 7t/2, ratio de cylindree 0,25. 

a. Evolution temporelle de la pression dans le tube a gaz. 

b. Evolution temporelle des debits massiques aux extremites. 

Figure 7. Piston phase shift 0.125 7t/2 to 0.8JT/2, with swept- 
volume ratio 0.25. a. Pressure temporal evolution inside tube 
of HPT refrigerator, b. Mass flow temporal evolution at hot 
and cold extremities of HPT refrigerator. 



Figure 5. Ratios de cylindree 0,25 ; 0,5 ; 0,75 ; phase 0,5K/2. 

a. Evolution temporelle de la pression dans le tube a gaz. 

b. Evolution temporelle des debits massiques aux extremites. 

c. Evolution temporelle du facteur de production frigorifique. 

Figure 5. Swept-volume ratios 0.25, 0.5, 0.75 ; phase shift 
0.5 71/2. a. Pressure temporal evolution inside tube of HPT 
refrigerator, b. Mass flow temporal evolution at hot and 
cold extremities of HPT refrigerator, c. Refrigerating factor 
temporal evolution of HPT refrigerator. 


Qf/(2Pi.fPV) Tf/Tc=0.8 



♦ Qf Vcpr=0 5 
o W Vcpr=0.5 
a Qf Vcpr=0.25 
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• Qf Vcpr=0.75 
o W Vcpr=0.75 


Figure 6. Evolution du facteur de production frigorifique et 
mecanique en fonction du dephasage des pistons du HPT, 
parametre en fonction du ratio des cylindrees. 

Figure 6. Refrigerating and mechanical factors evolutions for 
a HPT refrigerator depending on phase shift of pistons, with 
swept-volume ratio as parameter. 


Les performances presentees dans cet. article ne sont 
a considerer que de maniere relative ; il s’agit de valeurs 
maximales theoriques. obtenues avec un regenerateur 
d'efficacite E constante et unitaire. La modelisation 
du regenerateur est certainement une des taches les 
plus ardues qui s'offre au thermicien, sourtout avec 
les conditions de fonctionnement rencontrees dans les 
cryorefrigerateurs : temperatures et debits d'entree 
variables, compressibilite du fluide traversant la matrice 
poreuse. variation importante des proprietes physiques 
de celle-ci a basso temperature, etc. Pour Pavenir, il 
conviendrait done d'aborder ce probleme delicat. 

Malgre tout, la presente modelisation conjointe 
des deux types de refroidisseurs OPT et HPT fait 
ressortir, pour le second, deux parametres importants 
supplementaires, qui sont : 

- le rapport des cylindrees compresseur-dephaseur 

1 rep 

tip 

- la valeur de la phase entre les mouvements des 
pistons phas. 

L' influence antagoniste de ces deux facteurs impose 
de rechercher un reglage optimum du dephasage en 
fonction du rapport de cylindree choisi. Il semble 
preferable, pour la geometrie etudiee ici. de s’oriente-r 
vers une cylindree du module dephaseur netteinent 
inferieure a celle du compresseur (rapport 0,25), avec 
une phase voisine de 60 degres d’angle, cette valeur 
devant. etre diminuee pour un ratio superieur. La 
possibilite d’agir sur deux parametres, dephasage et 
cylindree, semble beaucoup plus souple au niveau des 
reglages de phase pression-debits que dans le cas du 
tube a orifice. 

Les performances theoriques calculees pour le HPT 
sont a peine inferieures a celles d un OPT. puisque Ton 
obtient un facteur de production frigorifique maximum 
de 0,03 contre 0.035, avec une efficacit.e frigorifique 

theorique cependant plus faible de 1.9 contre 

I, 

11,7 (AT r = 0,1) (figures 3d et. 6 ). Cependant, dans 
l’hypothese d'une recuperation du travail de detente 
du gaz dans le cylindre secondaire, on a calcule que 
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l'efficacite frigorifique peut etre ainelioree jusqu'a 5.56. 
Dans le cas d'un OPT. toute I'energie est degradee en 
pure perte dans rorifice. L'efficacite d : un OPT chute 
d'ailleurs tres vite si on augmente le facteur Ak r (sa 
valeur est de 0.5 pour Akr = 1). Celle du HPT varie 
assez peu pour les valeurs de dephasages autour du 
reglage optimal. 

En raison de relimination du reservoir, le HPT pour- 
rait offrir un encoinbrement moindre que le refroidis- 
seur classique OPT. des lors que des progres significatifs 
auront ete effect, ues sur les generateurs d onde do pres- 
sion. par exeinple en adaptant des microactionneurs 
piezo-elcctriques on autres. Ces avancces pcrmettraient 
d’envisager sa miniaturisation et son adoption defini¬ 
tive dans le domain© du refroidissemeut de coniposants 
electroniques ou d'appareils embarques. 
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Abrigded English Version 

Parametrical study for performances comparison between classic orifice pulse tube 
and hybrid Stirling pulse tube refrigerators 


This paper presents a parametrical theoretical 
comparative study of orifice pulse tube (OPT) and 
Stirling hybrid pulse tube (HPT) refrigerators. The 
HPT cryocooler uses a working gas under pressure 
enclosed in a variable volume. In the HPT refrigerator, 
two pistons, whose movements are generally non-phased, 
are disposed on each part, of the regenerator and 
pulse tube assembly. These piston movements allow 
gas cycle realisation and frigorific production. The HPT 
cryocooler is very similar to the Stirling refrigerator 
but with an expansion working piston at. ambient 
temperature. 

The internal dynamics of pulse tube refrigerator can 
be explained very succinctly in the following manner. 
When the compression piston moves to reduce the 


chamber volume, the gas is adiabatically compressed, 
and then, the compressor heat exchanger cools it 
to room temperature. One part of the compressed 
gas crosses the regenerator before entering the tube, 
tlie other part, goes out the tube to penetrate into 
the reservoir. A hot, heat exchanger, located between 
tlie tube end and the reservoir keeps the gas at 
ambient temperature and evacuates approximately the 
heat extracted in the cold heat exchanger. When the 
compressor piston expands, the expansion cools down 
the gas in the tube to a temperature lower than at 
the outlet, of the regenerator. Further in the expansion, 
the pressure becomes lower than the pressure of the 
reservoir: consequently, the gas flows out. of the reservoir 
and pushes the cold gas into the regenerator in the 
opposite direction. 
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Etude parametrique comparee des performances des refrigerateurs a tube a gaz pulse classique et hybride 


The frigorific production at the tube cold end 
is explained by axial enthalpy streaming [1-3], This 
stream, which has a temporal average equal to zero, 
is given by (1). It is maximal, in case of periodical 
variations of pressure and velocity, when this two 
variables are in phase. Figure 1 gives an idea of 
enthalpy distribution along various refrigerators. We can 
notice that, for an established periodical mode, average 
velocity of fluid is zero and regenerator gas temperature 
stays constant, if its thermal capacity is negligible 
compared to regenerator solid matrix capacity. So the 
enthalpy streaming is equal to zero in the regenerator. 

The comparison of the OPT and the HPT cryocoolers 
is based on results of a theoretical modelling. A system 
of dimensionless equat ions, based on mass conservation 
(4), ideal gas law (5) and energy conservation (C) is 
written for seven separate control volumes defined in 
the HPT by figure 2b. 

The two pistons movements (7) are supposed to be 
perfectly sinusoidal with adjustable phase shift: phas. 

In the regenerator, we make the hypothesis of li¬ 
near variations of pressure (8) and of solid matrix 
and gas temperature (9). We take into account the 
regenerator efficiency E to calculate regenerator gas 
outlet temperature. Mass flow rate through the regene¬ 
rator is represented by the following formula: 

rhdi = K 0 {Pc(t) - see also (19) 

Alternate heat transfer occurs between the fluid and 
the walls. Previously [15]. we have presented a formula¬ 
tion for this type of transfer, using compound numbers 
to take into account the phase between pressure, gas 
temperature, wall temperature and parietal heat flow. 
Real part of heat flow is: 

S f [ Qv dx = - ‘^1 Jfil L. (/> _ P) cos 3 
Jo VPe 7 - 1 

For mass flow rate calculations, we use the Poiseuille 
law for capillary tube, introducing pressure drop, length 
and diameter. The reservoir orifice is characterised by 
an aperture factor K^ T (13) and (18). So: 

m r = -A r (Prcp(t) - Pi{t)) 

Table I summarises the variables of the complete 
equations system: the total of balances gives seven 
differential equations for variables: 

Pc{t), Pt(t), Prcp(t), Tcp(t), Tt(O.f), T x (L.t), T rcp (t) 

Once these values are calculated, integration of seven 
other equations gives mass flow rate m c . m se* rh, cp , rhdi, 
m r . mi, rh, another equation for refrigerating power 
SQf and one for the work exchanged by the compressor 
piston: 

SIT = — +nic P T * 
dt 


Complete differential equations system is numerically 
solved, step-by-step, according to mass flow rate sign 
(inlet > 0. outlet < 0) and temperature choice (15). 
Table II shows the necessary eight successive phases for 
a complete cycle of working fluid. 

Written in dimensionless form, the equations empha¬ 
sise some important groups common for the OPT and 
the HPT and two additional ones typical for the HPT 
crvocooler: phase shift phas and second compressor vo¬ 
lume ratio V rC p/Vt. Others are: refrigeration production 
factor (16), volumes ratios: 


Pep I cl Pec Pc Pse 

Pt Pt Pt Pt Pt 


thermal factors (20) and flow of gas factors Akr and 
A'u 0 (18) and (19). First for the OPT. in figure 3, wo 
show that A'k, augmentation decreases the frigorific 
production, mainlv because the extreme pressures ratio 
p 

p“* x decreases. Pressure phase is also reduced and 

the two mass flow's of the tube extremities increase: in 
fact, their phases are almost constant. We also normally 
record that the best configuration occurs when pressure 
and flow rates are in phase (with small AW)- When 
AT, becomes too important, heat transfers of the tube 
extremities are inverted and the tube is transformed 
into an acoustical energy generator. 


For the HPT refrigerator in figures 4 and 5. we 

V rr , 

show an opposite influence of the two parameters 

Pt 

and phas, proving it is necessary, in order to obtain 
the greatest frigorific production, to research adequate 
phase shift for each working conditions. Generally, a 
smaller secondary swept-vohime is preferable for an 
efficient crvocooler (see figure 6). In contrast to the 
OPT. figure 7 proves the greater facility of the HPT 
crvocooler concerning amplitude and phase setting for 
rh c . The phas parameter lias a small influence on 
the pressure phase and the mass flow' rate rh s ,.. The 


theoretical maximal performance coefficient 


_SQ f 

P*V t 


of 


(HPT. = 0.25. phas 60 c ) is 0.03. a little less than 
Pt 


(OPT, A'kr = 0.1)) with 0.035. Its efficiency ^ is 1.9 

rr 


(5.56 if w r e count a possible recovery of expansion work) 
compared to 11.7. This efficiency falls rapidly with AW 
(value 0.5 for AW = 1), a parameter which is difficult to 
set precisely in practice. 


With further potential progress in reduced scale 
technology for pressure generators, for example by using 
piezo-electrics components and without a reservoir, the 
HPT crvocooler can have a lesser dimension than the 
OPT one for equivalent performances. We consider 
that in the near future the Stirling hybrid pulse tube 
crvocooler has good chances of being developed for 
miniature cryocoolers and spatial refrigerators. 
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